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In this paper we study thick-shell braneworld models in the presence of a Gauss-Bonnet term. 
We discuss the peculiarities of the attainment of the thin-shell limit in this case and compare them 
with the same situation in Einstein gravity. We describe the two simplest families of thick-brane 
models (parametrized by the shell thickness) one can think of. In the thin-shell limit, one family 
is characterized by the constancy of its internal density profile (a simple structure for the matter 
sector) and the other by the constancy of its internal curvature scalar (a simple structure for the 
geometric sector). We find that these two families are actually equivalent in Einstein gravity and 
that the presence of the Gauss-Bonnet term breaks this equivalence. In the second case, a shell will 
always keep some non-trivial internal structure, either on the matter or on the geometric sectors, 
even in the thin-shell limit. 
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I. INTRODUCTION 

In the light of the Randall-Sundrum (RS) braneworld paradigm pj, the cosmological evolution of the universe 
would acquire important modifications at early times with respect to the standard lore provided by the Friedmann 
equation of General Relativity 0, 0, (an exhaustive study of the phase space of the new cosmological models can 
be found in The new Friedmann-like equation has now contributions that are quadratic in the density, dark- 
radiation contributions and possibly even other types of dark contributions (coming for example from a fundamental 
electromagnetic field in the bulk [(| 0, @ ) . 

This same paradigm, with the characteristic presence of extra dimensions, naturally suggests that we should look 
for further modifications to the generalized Friedmann equation caused by the existence of a Gauss-Bonnet (GB) term 
in the field equations. In five dimensions, the most general geometric field equations of the form Gab — 0, with Gab 
a symmetric and conserved tensor constructed by using the metric and its first and second derivatives only, is a linear 
combination of the metric tensor itself, the Einstein tensor and the Lanczos tensor (or GB term) 0, 0] . From this 
point of view, the value of the three coefficients of this combination, basically, the cosmological constant, Newton's 
constant and the a constant (to be defined later) respectively, are parameters to be determined experimentally 1 . The 
GB term has particular relevance when considering the string-inspired nature of the RS proposal. In string theory, 
the GB term appears as the first higher curvature correction to Einstein gravity pd . 112. 1131 . Il4l| . Possible consequences 
of the GB interactions on cosmological inflation have already been considered in . Other analysis concerning 

different aspects of GB corrections to braneworld physics can be found in [Hl-lil (this is not intended to be an 
exhaustive list). 

In trying to generalize the thin-shell cosmological models to incorporate the effects of a GB term it has been some 
controversy. First, it was claimed that the GB equations were not well defined in the distributional sense required 
by the thin-shell formalism |17| . It would be necessary to study the microphysics to solve the ambiguities that would 
arise. However, later, in references ^^"Eli it was assumed that the particular structure of the GB Lagrangian (linear 
in second order derivatives) posed no problems in the distributional limit (at least in models with Z2 symmetry). 
Remarkably, the results obtained in those paper were not consistent with each other. The generalized Friedmann 
equation found in papers j^-J^I (quadratic on the brane density pi>) is different from that found in papers 
(with a complicated dependence on pj, coming from obtaining the real root of a cubic equation). More recently, there 
has been an argument |37ll38l| based on the form of the surface term in GB (see 39]) that strongly support this later 
complicated pb dependence (see also |4fj|). 

In this paper, we will study the cosmological behaviour of shells (or branes) that are thin but still of a finite thickness 
T. In this way we want to shed some light on how the zero thickness is attained in the presence of GB interactions. 
This limit has been studied for Einstein gravity in [4lJ. Thick shells in the context of GB interactions have been 
already studied in 42] and 43], but with a focus on different aspects than those in this paper. The conclusion of our 
analysis here is twofold. On the one side, our results show that the generalized Friedmann equation in ^^]-[24| can 



We are thinking on the standard RS scenario, in which there is only one additional extra dimension; for 7 or more extra dimensions one 
will have to consider more general Lovelock terms [j]. 



be found by using a completely general procedure , in which the energy density of the brane in the thin-limit is related 
to the averaged density. Instead, the equation in pRj— [27| can only be found for specific geometric configurations and 
with a procedure in which the energy density of the brane in the thin-limit comes from the value of the boundary 
density in the thick-brane model. On the other side, we argue that the information lost when treating a real thin shell 
as infinitely thin is in a sense larger in Einstcin-Gauss-Bonnet gravity than in the analogous situation in standard 
General Relativity. 

Let us explain further this last point. From a physical point of view, in the process of passing from the notion 
of function to that of distribution one loses information. Many different series of functions define the same limiting 
distribution. For example, the series 

f for \y\ >T/2, (0 for \y\ > T/2 , 

Mv) = { 9T(y) = { 2 (l) 

{ i for \y\<T/2; { for \y\<T/2, 

define the same limiting Dirac's delta distribution. The distribution only takes into account the total conserved area 
delimited by the series of functions. The gravitational field equations relate geometry with matter content. If we take 
the matter content to have some distributional character, the geometry will acquire also a distributional character. 
When analyzing the thin- limit of branes in Einstein gravity, by constructing series or families of solutions parametrized 
by their thickness, we observe that the blowing up parts of the series of functions that describe the density-of-matter 
profile transfer directly to the same kind of blowing up parts in the description of its associated geometry. Very simple 
density profiles [like the previous function friy)] are associated with very simple geometric profiles, and vice- versa. 
However, when considering Einstein-Gauss-Bonnet gravity this does not happen. The blowing up parts of the series 
describing the matter density and the geometry are inequivalent. A simple density profile does not correspond to a 
very simple geometric profile and vice- versa; on the contrary, we observe that they have some sort of complementary 
behaviour. This result leads us to argue that the distributional description of the cosmological evolution of a brane in 
Einstein-Gauss-Bonnet gravity is hiding important aspects of the microphysics, not present when dealing with pure 
Einstein gravity. Also, we find that for simple models of the geometry, one can make compatible the two seemingly 
distinct generalized Friedmann equations found in the literature. 

The paper is organized as follows. In the next Section we build thick-shell models for static branes (direct gen- 
eralization of the RS model ) ■ This provides a simple situation in which the main ideas of the paper can already 
be seen. For clarity, we will separate the Einstein and Einstein-Gauss-Bonnet cases. Section ITTT1 will deal with the 
dynamical cosmological case. Finally, we will make a brief summary of the results found in Section IIVI General 
formulae for the construction of the field equations are given in the Appendix. 

II. STATIC THICK SHELLS IN EINSTEIN AND GAUSS-BONNET 

A. Einstein 

To fix ideas and notation let us first describe the simple case in which we have a static thick brane in an anti-de 
Sitter (adS) bulk. We take an ansatz for the metric of the form 2 

ds 2 = e- 2A( -^r)^dxi J -dx u + dy 2 , (2) 

where r/^ is the four-dimensional Minkowski metric. Comparing with the formulas given in Appendix^[in particular 
equation IjAljl ] this means taking a(t,y) = n(t,y) = cxp(— 2A(y)), and b(t,y) = 1. The energy-momentum tensor has 
the form 

k§ T AB = Ruaub + Pl^ab + PrnAn B , (3) 

where ua = (— e 2A 7 0, 0), ua — (0,0, 1) and Hab = 9ab + uaub — nAns with gAB the 5-dimensional bulk metric. 
Here, k 2 denotes the bulk gravitational coupling constant and p, p^ and pr represent respectively the energy density, 
the longitudinal pressure and the transverse pressure, and are taken to depend only on y. The Einstein equations 



2 Five-, four-, and three-dimensional indices are written using upper-case Latin, Greek, and lower-case Latin letters respectively. 



Gab = —A-59AB + k 2 Tab with a negative cosmological constant, A5 = — 6/£ 2 , result on the following independent 
equations for the metric function A(y): 

3A"-6A' 2 = p-^, (4) 

6A' 2 =p T + ^, (5) 

Pl = -P- (6) 

For convenience, we will hide the k 2 dependence inside the matter magnitudes, p — «g|0truej etc. We also consider a 
Z2-symmetric geometry around y — 0. The brane extend in thickness from y — —T/2 to y = +T/2. Outside this 
region p = pt — 0, so we have a purely adS spacetime: A(y) = —y/£ + b for y G (—00, — T/2) and A(y) = y/£ + b for 
y G (T/2, +00) where b is an irrelevant constant. The junction conditions at y = —T/2,+T/2 [see Eqs. I|A6|) and 
(|A"7|) tell us that 

A(-T~/2)=A(-T+/2), A(T-/2) = A(T+ /2) , (7) 
A'(-T-/2) = A'(-T+/2), A'(T-/2) = A'(T + /2). (8) 

From here and using © we deduce that the transversal pressure is zero at the brane boundaries pt(—T/2) = 
Pt(T/2) = 0. Since we are imposing a ^-symmetry with y — as fixed point, hereafter we will only specify the value 
of the different functions in the interval (— T/2,0). 

The function A' is odd and therefore interpolates from A'(— T/2) = — 1/£ to A'(Q) = 0. If in addition we ask for the 
null-energy condition p + pr — 3A" > to be satisfied everywhere inside the brane, then pt has to be a negative and 
monotonically decreasing function from pt(—T/2) = to pt(0) = —6/£ 2 - This condition will turn to be fundamental 
in defining a thin-shell limit. 

By isolating A" from equations (0J and (|5j) we can relate the total bending of the geometry on passing through the 
brane with its total p + pt 

^=3A'\ T g /2 = J T '\p + p T )dy. (9) 

At this stage of generality, one can create different one-parameter families of thick-brane versions of the Randall- 
Sundrum thin brane geometry, by parameterizing each member of a given family by its thickness T. The only 
requirement needed to do this is that the value of the previous integral must be kept fixed independently of the 
thickness of the particular thick-brane geometry. Thus, each particular family can be seen as a regularization of 
Dirac's delta distribution. 

We can realize that, provided the condition p + pt > is satisfied, there exists a constant C, independent of the 
thickness T, such that pt < C ', that is, the profile for pt is bounded and will not blow up in the thin-shell limit. 
Therefore, in the limit in which the thickness of the branes goes to zero, T — > 0, the integral of px goes to zero with 
the thickness. (Strictly speaking, the thin-shell limit is reached when T/£ — * but throughout this paper we are going 
to maintain £ constant.) Instead, the profile of p has to develop arbitrarily large values in order to fulfil 

r T/2 



6 



£ T^O 



lim / p dy. (10) 



T/2 



In the thin-shell limit, we can think of Einstein's equations as providing a relation between the characteristics of the 
density profile and the shape of the internal geometry. A very complicated density profile will have associated a very 
complicated function A(y) . Physically we can argue that when a shell becomes very thin one does not care about its 
internal structure and, therefore, one tries to describe it in the most simple terms. But what it is exactly the meaning 
of simple? Here we will adopt two different definitions of simple: The first one is to consider that the internal density 
is distributed homogeneously throughout the shell when the shell becomes very thin. The second case is to consider 
that the profile for A' is such that it interpolates from A' (—T/2) = —l/l to A'(Q) = through a straight line, or what 
is the same, that the internal profile of A" is constant. Again, we require this for very thin shells. This geometric 
prescription is equivalent to asking for a constant internal scalar curvature, since R = 8^4" — 20A' 2 and for every thin 
shell the term A' 2 is negligible with respect to the constant A" term. Hereafter, we will use indistinctly the names 
straight interpolation or constant curvature for these models. In building arbitrarily thin braneworld models, one 
needs that the profiles for the internal density p and the internal A" acquire arbitrarily high values (they will become 
distributions in the limit of strictly zero thickness). In the first of the two simple models described, the simplicity 
applies to the blowing up parts of the matter content; in the second, the simplicity applies to the blowing up parts of 
the geometry. From the physical point of view advocated in the introduction, these simple profiles are those that do 
not involve losing information in the process of taking the limit of strictly zero thickness. 

Let us analyze both cases independently. 



1. Constant density profile 



Let us first define for convenience z = y/T as a scale invariant coordinate inside the brane. Then, mathematically, 
the idea that the density profile, which we will assume to be analytic inside the brane for simplicity, becomes constant 
in the thin-shell limit can be expressed as follows: 

p(z) = ^/3 n (T)z 2 ", (11) 

TO 

where 

Urn T l3 n {T) -> 0, Vn ^ 0; hm T /? (T) -> p b : constant . (12) 
For these density profiles, the Einstein equations in the thin-shell limit tell us that 

3A" =p (T)-^+6A' 2 . (13) 

From here we get the profile for A': 



Notice that this expression only makes sense for Pq(T) > 6/£ 2 , but this is just the regime we are interested in. We 
have to impose now the boundary condition A'(T/2) = l/l to the previous expression (|14[) 



1 f3 (T) 1 / //3o(T) 1 \ 

€ = V - 6 ^ taJ1 ^V - 6 ^ T J • (15) 

In this manner, we have implicitly determined the form of the function /3q(T). In the limit in which T — ► with 
Tf3 {T) — > /%, we find the following relation 

^=P6- (16) 
This condition is just what we expected from the average condition i|l(J|) . 



2. Straight interpolation 

In this case, the mathematical idea that in thin-shell limit the profile for A' corresponds to a straight interpolation, 
can be formulated as 

A"(z)=J2^(T)z 2n , (17) 

n 

with 

limT 7n (T)^0, Vn^O; lim T 7o (T) -» ^ (18) 

For these geometries we find that the associated profiles for px and p in the thin-shell limit have the following form 

PT = -Jr(l-4z 2 )+c^(T,z), (19) 

/°=^ + |( 1 - 4 ^)+^(T^) ) (20) 

where here and along this paper u) n (T, z) will denote functions that vanish in the limit T — > 0. Now, from this density 
profile we can see that 



fT/2 g 

hm / pdy= 7 , (21) 

-T/2 l 



T->0 



as we expected. Moreover, we can see that the boundary value of the density satisfies Tp\ T / 2 ~^ 6/£ in the thin shell 
limit, which is the same condition satisfied by the averaged density, T(p) — > 6/£. 

An additional interesting observation for what follows is the following: The sets of profiles that yield constant density 
in the thin-shell limit and straight interpolation for the geometric profile H17|) coincide. Therefore, in the thin- 
shell limit one can assume at the same time a constant internal structure for the density and a straight-interpolation 
for the geometry. 

B. Einstein-Gauss-Bonnet 

Let us move now to the analysis of the same ideas in the presence of the Gauss-Bonnet term. The Einstein-Gauss- 
Bonnet Lagrangian is 

S=^J dx 5 ^[R-2A 5 + aL GB ], (22) 

with 

Lgb = R 2 — ^R ab Rab + R ABCD Rabcd ■ (23) 
Now, the field equations deduced from this Lagrangian are 

Gab + olR ab = -A 5 g AB + k\T ab , (24) 

where Hab is the Lanczos tensor [T(l |: 

Hab = 2RacdeRb DE ~~ ^RacbdR CD — ^RacRb + 2RRab — -^QabLgb ■ (25) 

For the ansatz J5J we obtain [compare with equations l|A18IA20(l in Appendix E] 

3A"(1 - AaA 12 ) - &A' 2 {1 - 2aA' 2 ) = p-%, (26) 

6A' 2 (l-2aA' 2 )=p T + ^, (27) 

PL = -P- (28) 

The junction conditions for the geometry are the same as before JSJ, implying again the vanishing of the transversal 
pressure at the boundaries, px = 0. 

In the outside region the solution is a pure adS spacetime but with a modified length scale 
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Now, isolating A" from l|2tj[l and (|27|l we can relate the total bending of the geometry on passing through the brane 
with the integral of p + pt 

I f 1 - = $ A ' - 4aA' 3 )| T/ T 2 /2 = [ m (p + PT ) dy. (30) 

-T/2 



3pJ y n - T ' 2 



Again, if the condition p + px > is fulfilled throughout the brane we will have that in the thin shell limit 



6 / 4 a 



T/2 

i ( 1 ~ oi; ) =J, im / P d v- ( 31 ) 

T/2 



3 p J T^O 



At this point we can pursue this analysis in the two simple cases of constant density profile and straight interpolation. 



1. Constant density profile 



Following the same steps as before for a constant density profile (|f f () - (|f 2f> . the equation that one has to solve in the 
thin-shell limit is 



3A"(1 - AaA' 2 ) = A,(T) - |r + %A' 2 (l - 2aA' 2 ) . 
Introducing the notation B = A' we reduce this equation to the following integral 



4a 



The result of performing the integration is 



(AaB 2 -VjdB 



o B 4 - ±B 2 

2a 2a V 



: tan 



B 



where R± are 



Again, by imposing the boundary condition 

1 



(21 L 

6 e 2 



■. tanh 1 



(32) 



(33) 



B 



1± 4 1 + 8q 



T _ 1 
2 ~ 2 



■ tan 



Po(T) 1 
6 £ 2 



: tanh 



(34) 



(35) 



(36) 



we find the appropriate form for (3q (T) . With a lengthy but straightforward calculation we can check that in the limit 
T -> 0, /3 (T) -► oo, we have 



(i 



i 



4 a 



3 P 



(37) 



in agreement with condition 1)31(1. 

Using this same asymptotic expansion, we can see that, in the thin-shell limit, the profile for A'(y) satisfies 

Recursively, one can create a Taylor expansion for A'(y). The first two terms are 

A'{y) = ±/3o(T) y + ^/3o(T) 3 y 3 + 0(y 5 ) = l -T (3 q {T) z + ^T 3 /? (T) 3 z 3 + 0(z 5 ). 
By differentiating this expression we find 

A"(y) = ^/3o(T) + ^T 2 /? (T) 3 z 2 + G(z 4 ). 



(38) 



(39) 



(40) 



Now, contrarily to what happens in Einstein theory, this profile does not correspond to the set considered in the 
straight interpolation before (see Fig. 1). By looking at 1)171) we can identify 



7o (T) = ±A,(T), 7i (T)ee^T 2 /? (T) 3 . 



(41) 



Then, we can see that 



f 



4 a 



3^ 



32a 



limT 7 o(T) = 7 limT 7l (T) = ^- 1--^ ^0 



T^0 



4 a 



3 £2 



(42) 



The coefficients 7n do not satisfy the conditions in 1)18)1 . Therefore, unlike the energy density, the scalar of curvature 
does not have a constant profile. 



P(y) 




FIG. 1: Family of constant density profiles with decreasing thickness and associated geometric profile for A' and A". 



2. Straight interpolation 
As in the Einstein case, the straight interpolation profile for A" corresponds to 

A"(z)=J2^(T)z 2n , (43) 

n 

with 

HmT 7n (T)-0, Vn^O; UmT 7o (T) - |. (44) 

From here we can deduce the associated profiles for j?t and p by substituting on l(26|) and 1)27(1 . 
In the limit T — > 0, the dominant part in the density profile is 

p = 3 7o (T)(l - 4a 7o (T) 2 T 2 z 2 ). (45) 

Identifying 

/3o(r) = 3 7 o(T), /? 1 (T) = -12aT 2 7o (T) 3 , (46) 

we find that 

lim T (3 (T) = 1, lim T/3UT) ^ . (47) 

Therefore, even in the thin shell limit a straight interpolation in the geometry does not correspond to a constant 
density profile (see Fig. 2). In the presence of a Gauss-Bonnet term it is not compatible to ascribe to have a simple 
description for the interior density profile and for the geometric warp factor at the same time. In the limit of strictly 
zero thickness (distributional limit) one will unavoidably lose some information on the combined matter-geometry 
system. 

To finish this section let us make an additional observation. From expressions (|43|l and (|44p. we can see that 

A" = ~A'\ T/2 + u(T, z) with Hm2V(T, z) = . (48) 



Using this property in (|26I27|I 



p + p T = 3A"(1 - AaA' 2 ) = ^U'(T/2)(1 - 4aA' 2 ) (49) 




FIG. 2: Plot of the straight interpolation profile for the geometric factor A' and its associated density profile. 



and evaluating at y = T/2 we find 

Tp\ T/2 = 6A'(l-4aA' 2 )\ T/2 . (50) 
We can see that contrarily to what happens in the Einstein case, this condition is different from the averaged condition 



T(p) = lim / ' pdy = (3A' - AaA' 3 )\l^ /2 = 6A'(l - ^aA' 2 ) 

T ^°Jt/2 ' 



(51) 

T/2 



Therefore, the averaged density and the boundary density are different, and this is independent of the brane thickness. 
For this simple model, in thin-shell limit one can define two different internal density parameters characterizing the 
thin brane. One represents the total averaged internal density and can be defined as 

p av = lim T(p). (52) 

The other represents an internal density parameter calculated by extrapolating to the interior the value of the density 
on the boundary. This density can be defined as 

p hv = fim Tp\ T/2 . (53) 

The junction conditions for a thin shell given in 0], corresponds to the averaged condition lj31|) or (|51|l and therefore 
relate the total bending of the geometry in passing through the brane with its total averaged density. Instead, the 
particular condition analyzed for the boundary value of p\t/2 <|5U|I yields in the thin-shell limit the junction condition 
in |25j . This condition is only considering information about the boundary value of the density and not about its 
average. 

As a summary, what these analysis suggest is that in the presence of the Gauss-Bonnet term we cannot forget the 
interior structure of the brane, by modelling it by a simple model, even in the thin-shell limit. We will see again this 
feature in the next section on the cosmological dynamics of thick shells. 



III. DYNAMIC THICK SHELLS IN EINSTEIN AND GAUSS-BONNET 



We are now going to study cosmological thick branes. To that end we will use the class of spacetime metrics 
given in (|A1(I . which contain a Friedmann- Robertson- Walker (FRW) cosmological model in every hypersurface {y — 
const.}, with a matter content described by an energy-momentum tensor of the form Pjl. We consider the additional 



assumption of a static fifth dimension: 6 = 0. We can rescale the coordinate y in such a way that 6=1. Then, the 
line element (IA1(I becomes 



ds 2 = -n 2 (t,y)dt 2 + a 2 (t,y)hijdx l dx j + dy 2 



(54) 



In Appendix 1X1 we show that the {ty}— component of the Einstein-Gauss-Bonnet field equations, for the case with a 
well-defined limit in Einstein gravity, leads to the equation i|A13(l . In our case it implies the following relation: 



n{t,y) = £{t)a{t,y) 



(55) 



In this situation the rest of field equations can be written in the form given in (|A18IA19IA20(I . In our case they 
become 3 



a 4 $ + 2a$ 2 



(2 



(a 4 )V, 



(56) 



n a 
n a' 



$ + 2a$ 2 



$ + 2a$ 2 



2aa! a 2 pL 



(57) 



where now $ is given by 



a 4 $ + 2a$ 2 



k 



r(a 4 )'PT- 



a' 2 , k a' 2 

— TT I _ 



(58) 



(59) 



where we have identified the first term with the square of the Hubble function associated with each y — const, slide, 

H(t,y) = —. (60) 



With the assumption 6 = 0, the field equation l|57() leads to a conservation equation for matter of the same form as 
in the FRW models [see Ea. (|ATo)l ]: 



P = -3-0 + Pi) 
a 



(61) 



In the same way we did in the static scenario, we consider here the situation in which there is a Z 2 symmetry and 
a fixed proper thickness T for the brane. Then one has to solve separately the equations for the bulk (jy\ > T/2) and 
the equations for the thick brane (\y\ < T/2). The first step has already been done and the result is |2q : 



1 



M 



$ + 2a<P z + — = — 



for \y\ > I 



T 
2 " 



(62) 



where M is a constant that can be identified with the mass of a black hole present in the bulk. Once the solution 
inside the thick brane has been found one has to impose the junction conditions l|A6IA7|) at y = ±T/2. 

The first thing we can deduce from the junction conditions is that the quantity $ is continuous across the two 
boundaries y — ±T/2 . But in general, its transversal derivative, 4>', will be discontinuous. Then, using equation 158(1 
it follows that the transversal pressure has to be zero on the boundary, pr(t, ±T/2) = 0. At the same time, from l(58|) 
we deduce that we must always have: 



a 4 { $ + 2a$ 2 + - 



= M. 



(63) 



y=±T/2 



3 The coupling constant a used here is one half the one used in (25(1 . 



On the other hand, using again the relation (|A13|> we find that 

H' = --H => $' = -2-(W— V (64) 
a a \ a J 

and then, expanding 1)56(1 . we get 

(l + 4a<J>)^ = $-|-ip. (65) 

In the limit T — > 0, the profiles for the density p and for a" blow up, therefore, these dominant terms in expres- 
sion 165(1 have to be equated. This results in the following equation 

(l + 4a<J>)(jQ =-\p. (66) 
In what follows we consider the analysis of the Einstein and Einstein-Gauss-Bonnet theories separately. 



A. Einstein 



In Einstein gravity it is not difficult to write down an equation describing the dynamics of every layer in the interior 
of a thick shell. To that end we will take a = in the equations above. By integrating 1(56(1 over the interval (—T/2, y*) 
and using I(63I66|) we arrive at 



1 



M 
~7a 



l 

6a 1 



(« 4 )Wy = 4 

-T/2 6Q 



pdy 



M 

~7a 



l 

6a 1 



T/2 



(a 4 )' pdy. (67) 



A particular layer of matter inside the shell, located at y = y*, can be seen as separating an internal spacetime from 
a piece of external spacetime. From the previous equation we can see that the cosmological evolution of each layer 
y = in the thick shell depends on the balance between the integrated density beyond the layer (external spacetime) 
and a weighted contribution of the integrated density in the internal spacetime. Therefore, the dynamics of each shell 
layer will be influenced by the particular characteristics of the internal density profile inside the shell. However, by 
looking at this same equation we can see that the dynamics of the boundary layer = — T/2 is only influenced by 
the total integrated density throughout the shell: 



k 



1 



M 

~7a 



-T/2 




(68) 



This is the modified Friedmann equation for the cosmological evolution of the brane . 

In the same manner as we proceeded with static shells, let us analyze the case in which the density profile tends to 
a (time-dependent) constant in the thin-shell limit: 



p = 0o(T,t)+U3(T,t,z), 



lim T0o(T,t) 



lim Tu>3(T, t, z) 



= 0. 



(69) 



When a = 0, equation ((66(1 tells us that if the density profile depends only on t in the thin-shell limit, then, in this 
same limit, the blowing up part of the geometry (a! /a)' is also constant through the brane interior, describing what we 
called before a straight interpolation. A simple density profile amounts to a simple and equivalent geometric profile 
and vice- versa. In this same case but including the Gauss-Bonnet term, a / 0, the geometrical factor (a' /a)' will 
exhibit a non trivial profile in y, even in the thin-shell limit. We will see this fact in more detail in the next subsection. 
Now, in the case in which p depends only on time, equation 1(67(1 reads 



tt2 k 1 
H + — + 72 

From equation 1(66(1 we deduce that 
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(71) 



and integrating we obtain 



i(t,y) = a (t)exp ( --p av (t)Tz 2 ) , 



(72) 



(remember that z = y/T). Therefore, in the lowest order in T we have an equation for the internal geometry of the 
following form 



, k 1 loo lo, 9x M 1 
which is exactly the standard braneworld generalized Friedmann equation . 



(73) 



B. Einstein-Gauss-Bonnet 



In the general Einstein-Gauss-Bonnet case, equation can be written as 
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Then, integrating between —T/2 and T/2 yields 
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The boundary equation l|63[) can be written as 
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This is a quadratic equation for (a'/a) 2 |r/ 2 with solutions 
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From these two roots we will take only the minus sign as it is the only one with a well defined limit when a tends 
to zero. Now, by squaring l(75|l and substituting the previous root we arrive to a cubic equation for H 2 + k/a 2 first 
found in [l^ . This cubic equation has a real root that can be expressed as |44| 
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In addition to this equation we need the conservation equation 

P = 3H{p + p L ) , 



(78) 



(79) 



(80) 



which is valid for each section y = y*, and in particular, for the boundary, y = T/2. This equation can be averaged 
to give 



or written in another way 



T(p) = 3(HT(p + PL )}= 3H\ T/2 (T(p) + T(p L }) + 0(T), 



P&v = 3H\ T / 2 (p av + p avL ) 
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(82) 



This happens because 



ff(t,y)-> H (t,y ) + O(T) 



(83) 



for whatever yo € [— T/2, T/2], which we have taken as yo = r/2 for convenience. 

Let us analyze now the simple case of a constant density profile (|69p. For consistency with the T — > case, we 
know that 



a(t,y) = a (t)[l + T&(t,z)] + O(T 2 ). 

and therefore, from l(55jl 

n(t,y)=t(t) [ao(l + Ta(t,z))Y . 
Then, in the same limit equation l|66l) reads 
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(Here the subscript , z denotes differentiation with respect to z.) A necessary condition to have a straight interpolation 
for the geometry is that a(t, z) — b(t)Z(z). To check whether or not a simple density profile corresponds to a simple 
geometrical profile we can therefore try to solve the previous equation by separation of variables. It is not difficult to 
see that in order to find a solution with a well defined Einstein limit we need that 



b(t) = fi, n- l 0o{T,t)T = ^ _1 /3 (T)T = p w : constant, 
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(87) 



where [i is a constant that can be absorbed into the function Z(z), so we will take it to be fi = 1. In this way we will 
recover the anti-de Sitter and de Sitter solutions for the brane (depending on the sign of the effective four-dimensional 
cosmological constant). To find the specific y profile we have to solve 
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This equation can be integrated to get 
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(89) 



For our proposes the specific solution of this cubic equation is not important. What we want to point out is that the 
solution does not correspond to a straight interpolation as it happened in the Einstein case. So, in general, simple 
solutions for the matter profile lead to non trivial profiles for the scalar curvature even in the thin-shell limit. 
Let us see now what happens when we take a simple model for the geometry, the straight interpolation model: 



a(t,z) = ao(t) — —b(t)z 2 T. 
From equation l|74(l we can deduce the density profile in this situation 
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As in the static case, even for very small thickness the density profile has now a non-trivial structure. We can observe 
that 
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The second relation and equation (|48|l coincide in the thin-shell limit. Therefore, evaluating 174(1 on y = T/2 we 
obtain 
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(93) 



Now, following the same steps that we followed previously but using this condition instead of i|75|) we arrive at a 
cosmological generalized Friedmann equation pfij different from that in [Tsj in its form and in the fact that it depends 
on the quantity associated with the boundary value of the energy density, pbv, instead of the value associated with the 
average of the energy density, p av . Remarkably, the cubic equation that results from combining the last equation l|9"5|) 
with the boundary condition (|76|l becomes in this case linear. That is, the coefficients of the terms quadratic and 
cubic in H 2 + k/a 2 vanish |2£j. The modified Friedmann equation found in this case is: 

, k 1 1 2 1 / r 8a 8aM \ 

H ~ + ^ = ^TW^W)^ Phv + ^{^ — - 1 ) ■ (94) 

In contrast with the modified Friedmann equation (|78|l . which was obtained by using a completely general procedure, 
in order to obtain this equation we had to use a procedure which required to consider an extra assumption, namely 
equation (|92|) . and hence, it will not work for profiles of the metric function a(t, z) that do not satisfy these requirements 
or equivalent ones. On the other hand, by looking at the developments here presented, we can conclude that the 
different results found in the literature for the dynamics of a distributional shell have their origin in the additional 
internal richness introduced in the brane by the presence of the GB term. 



IV. SUMMARY 



We have analyzed and compared how the thin shell limit of static and cosmological braneworld models is attained 
in Einstein and Einstein-Gauss-Bonnet gravitational theories. We have seen that the generalized Friedmann equation 
proposed in is always valid and relates the dynamical behaviour of the shell's boundary with its total internal 
density (obtained by integrating transversally the density profile). Instead, the generalized Friedmann equation 
proposed in |25| relates the dynamical behaviour of the shell's boundary with the boundary value of the density 
within the brane. This equation is not always valid, only for specific geometrical configurations. 

Einstein equations in these models transfer the blowing up contributions of the thin-shell internal density profile to 
the structure of the internal geometry in a faithful way. If we don't know the internal structure of the shell we can 
always model it in simple terms by assuming an (almost) constant density profile and an (almost) constant internal 
curvature. However, the GB term makes incompatible to have both magnitudes (almost) constant. If the density is 
(almost) constant, then the curvature is not, and vice-versa. Therefore, we can say that the particular structure of 
the Einstein-Gauss-Bonnet theory introduces important microphysical features on the matter-geometry configurations 
beyond those in Einstein gravity, that are hidden in the distributional limit. 
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APPENDIX A: THE 5D METRIC AND ITS ASSOCIATED GEOMETRICAL QUANTITIES AND FIELD 

EQUATIONS 

In this appendix we present the main geometrical quantities and field equations associated with the following 5D 
metric 

ds 2 = g A B<ix A dx B = -n 2 (t 1 y)dt 2 + a 2 (t 1 y)h ij {x k )dx i dx j + b 2 {t 1 y)dy 2 , (Al) 

where hij is the metric of the three-dimensional maximally symmetric surfaces {t — const., y = const.}, whose spatial 
curvature is parametrized by k = — 1,0, 1. A particular representation of is 

h lj dx l dx j = - j (dr 2 + r 2 dnj) , (A2) 

(1 + §r 2 ) 

being dfl 2 , the metric of the 2-sphere. The metric IjAlfl contains as particular cases the metrics used along this paper. 



The non-zero components of the Einstein tensor Gab corresponding to this line element are given by (Q = dtQ, 
Q' = d y Q): 
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Apart from the metric and the Einstein tensor, the held equations in Einstein-Gauss-Bonnet gravity (|24|1 contain a 
term quadratic in the curvature, namely Hab [see Eq. I)25[l]. The non-zero components of this tensor can be written 
as follows 
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(A5) 



In this paper we consider the situation in which a thick brane is embedded in the five-dimensional spacetime 
described by IjAlj) . whose boundaries are located at y = const, hypersurfaces. Let us consider the usual junction 
conditions at a hypersurface E Wc = {p G V5 \ y(p) = y c }, that is, the continuity of the induced metric, qab 



9ab — riAns and the extrinsic curvature, Kab = — ^{aIb)^ C n Di °f 



n(t,y+) = n(t,y c ) 
n'(t,y+) n'(t,y-) 



a{t,yt) = a(t,y c ) j 
a'(t,y+) a'(t,y-) 



b(t,y+) b(t, V o) ' b(t,y+) b{t, V c) ' 
Now, let us assume a matter content described by an energy-momentum tensor of the form 
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«5 Tab = Puaub + Vl^ab + pt"-a"-s , 



(A8) 



where 



u A = (-n(t,y), 0,0) , h AB = g AB + u A u B - n A n B , n A = (0, 0, b(t, y)) , 



(A9) 



where p , pi, , andpr denote, respectively, the energy density and the longitudinal and transverse pressures with respect 
to the observers u A . They are functions of t and y. The energy-momentum conservation equations, V A T AB = 0, 
reduce to the following two equations: 
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In this situation, the {t?/}-component of the field equations for the metric (jAip in Einstein-Gauss-Bonnet gravity 
[Eq. (|23}] has the form 
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If we discard the possibility 1 + 4a$ = by restricting ourselves to models with a well-defined limit in Einstein gravity 
(a — » 0), we have that the metric functions must satisfy the following relation 
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Using this consequence of the {iy}-component, we can rewrite the rest of components of G AB and H AB as 
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Then, the field equations Ij24(l for the metric IjAip are equivalent to equation l|A13jl and the following three equations 



a 4 [ $ + 2a$ 2 + p 



g(« 4 )V, 



(A18) 



6 d 



a 4 $ + 2a$ 2 + 



£ 2 



n' a 
n a' 



a 4 $ + 2a$ 2 



£ 2 



a 4 $ + 2a$ 2 



£ 2 



2aa'a A p L , (A19) 



« f <^> + 2a^ + - 



g(a 4 )'PT- 



(A20) 



Introducing (|A18|I and l|A20|) into i|A19|) we get the conservation equation (|A10|I . 
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